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ABSTRACT 


A  number  of  theoretical  calculations  of  the 
morion  of  a  spherical  underwater  explosion  bubble  oscil¬ 
lating  in  an  incompressible  homogeneous  unbounded  in- 
viscid  fluid  are  presented,  using  various  "models"  for 
the  behavior  of  the  bubble  interior.  The  purpose  of 
this  study  was  to  find  an  adequate  representation  for 
U9e  in  axisymmetric  calculations  of  underwater  nuclear 
explosion  effects  in  which  the  bubble  may  become  non- 
spherical  due  to  the  effects  of  gravity  and/or  nearby 
boundaries.  The  cases  considered  include  (1)  neglect¬ 
ing  the  bubble  atmosphere  altogether;  that  is,  trea  ' 
the  bubble  as  an  evacuated  cavity,  (2)  treating  tne 
bubble  interior  as  an  adiabatic  homogeneous  ideal  gas, 
but  ignoring  the  inertial  effects  of  the  gas,  (3)  in¬ 
cluding  the  gas  inertial  effects  in  the  previous  case, 

(4)  replacing  the  ideal  gas  behavior  in  case  2  with 
equat ion-of-state  data  for  real  steam,  and  (5)  revert¬ 
ing  to  an  ideal  gas  for  simplicity ,  dropping  the  assump¬ 
tion  of  bubble  homogeneity;  that  is,  the  gas  dynamics  of 
the  buoble  interior  were  followed  by  explicit  numerical 
integration  of  Euler 'b  equations  using  a  Lagrangian 
finite-difference  computer  code.  Except  for  case  1, 
the  results  were  found  to  be  virtually  identical  and  in 
good  agreement  tdth  with  experimental  measurements  and 
with  calculations  performed  by  other  investigators  which 
take  into  consideration  the  compressibility  of  the  water 
outside  the  bubble. 
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I .  INTRODUCTION 

Underwater  explosion  effects  have  been  the  subject  of 
extensive  study  for  many  years,  particularly  since  the  beginning 
of  World  War  II.  Most  of  the  early  research  concerned  the 
characterization  of  the  shock  wave  caused  by  the  explosion  find 
its  interaction  with  the  water  surface,  the  aea  bottom,  and 
(especially)  targets  such  as  the  hulls  of  ships  and  submarines. 

The  effects  of  the  shockwave  may  be  felt  at  very  lar^e  distances 
and  occur  on  a  very  short  time-scale  due  to  the  high  speed  of 
sound  in  water.  The  late-time  effects  of  the  explosion  (that  is, 
the  pulsation  and  migration  of  the  residual  "bubble"  of  gas, 
its  eruption  from  the  surface,  the  resulting  water  waves,  and 
similar  phenomena)  received  much  less  attention.  This  was,  of 
course,  only  natural,  since  the  military  importance  of  these 
phenomena  is  usually  relatively  marginal  for  conventional  under¬ 
water  ordnance. 

When  nuclear  (and  later,  thermonuclear)  explosives 
were  developed,  however,  the  available  energy  release  of  under¬ 
water  weapons  increased  by  several  orders  of  magnitude.  The 
secondary  late-time  mass-motion  effects  consequently  acquired 
new  significance,  and  therefore,  in  recent  years,  these  effects 
have  been  much  more  extensively  studied  than  previously.  One 
of  the  effects  of  underwater  nuclear  explosions  which  immediately 
comes  to  mind  is  that  of  the  dispersal  of  the  radioactive  nuclear 
deDris;  the  transport  of  this  residue  from  the  point  of  burst  to 
the  above-surface  environment  is  governed  by  tha  late-time  motion 
of  the  steam  bubble  produced  by  the  explosion.  Another  effect 
of  military  significance,  particularly  for  large  explosions,  is 
thw  generation  of  large  water  surface  waves  whioh  could  pose  a 
serious  threat  to  naval  units  or  could  even  inundate  a  nearby 
coastline . 


The  mass-motion  effects  of  an  underwater  explosion,  as 
compared  to  these  cf  tha  shockwave,  ar#  fairly  slow  an  'eiatant, 
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After  the  emission  of  the  shock,  the  hot  gases  formed  by  the 
explosion  expand;  if  the  explosion  is  shallow  enough,  this  gas 
bubble  may  erupt  from  the  water  surface  during  the  first  expan¬ 
sion.  Otherwise,  the  bubble  will  expand  to  a  maximum  size,  then 
contract  to  a  minimum  size,  re-exp*n<l,  and  continue  to  pulsate 
with  diminishing  amplitude.  At  the  same  time,  the  bubble  will 
migrate  upward  toward  the  surface  due  to  its  buoyancy.  For 
nuclear  explosions,  thlB  bubble  oscillation-migration  phase 
occurs  on  a  time  scale  of  several  seconds  as  compared  to  milli¬ 
seconds  for  shock  wave  effects,  and  therefore  the  motion  is 
virtually  always  subsonic.  The  eruption  of  the  bubble  from  the 
surface  will  hurl  large  masses  of  water  aloft,  releasing  the 
contained  fission  products,  generating  large  surface  waves  that 
thereafter  propagate  away  from  surface  zero,  and  setting  into 
motion  other  familiar  late-time  explosion  effects  (the  residual 
upwelling  along  the  explosion  axis,  the  turbulent  diffusion  of 
the  radioactive  surface  "pool",  and  so  on). 

Since  the  bubble  pulsation,  migration,  plume  eruption 
and  subsequent  events  are  fairly  slow-motion  phenomena,  the 
water  motion  may  be  adequately  treated  as  incompressible  flow. 
Even  so,  however,  the  problem  is  in  general  at  least  tvo-dimen- 
sior'l  (that  is,  axisymmetric)  and  time  dependent,  and  involves 
f’-ee  surfaces.  Therefore,  the  theoretical  treatment  of  the 
motion  is  quite  complicated,  and  consequently  most  of  the  infor¬ 
mation  available  today  has  been  collected  by  experimental  means. 
Recently,  however,  the  development  of  modern  high-speed  third 
generation  digital  computers  has  rendered  purely  theoretical 
calculations  of  axisymmetric  bubble  motion  feasible. 

The  MACYL6  hydrodynamic  code  (Pritchett,  1970a)  was 
designed  specifically  to  compute  the  water  motion  around  a 
pulsating,  migrati-g  explosion  bubble,  including  late-time 
plume  eruption  and  subsequent  phenomena,  by  brute-force  numerical 
integration  using  finite-difference  techniques  of  the  fundamental 
governing  equations  of  hydrodynamics  in  axisymmetric  geometry. 
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The  computer  program  has  been  used  successfully  in  the  past  to 
compute  the  late-time  flow  after  various  underwater  nuclear 
explosions,  and  agreement  with  experimental  results  (where 
available)  has  been  good  (see,  for  example,  Pritchett  and 
Pestaner ,  1969;  Pritchett,  1970b;  Pritchett ,  1971).  It  is  now 
being  put  to  use  to  study  the  formation  of  surface  waves  by 
deep  underwater  nuclear  explosions. 

"n  order  to  use  the  MACYL6  code  to  describe  an  under¬ 
water  explosion,  however,  the  boundary  conditions  must  be 
specified;  air  pressure  is,  of  course,  imposed  at  the  air-water 
interface  and  remains  constant  with  time.  At  the  bubble-water 
interface,  a  pressure  must  also  be  prescribed  in  some  realistic 
manner.  That  is,  the  MACYL6  code  computes  the  flow  in  the  water 
surrounding  the  bubble,  but  the  gas  dynamics  within  the  bubble 
itself  must  be  "modelled"  in  some  approximate  fashion  so  as  to 
supply  the  required  boundary  pressure.  In  the  calculations 
published  to  date,  the  bubble  pressure  was  assumed  to  be  a  func¬ 
tion  of  bubble  volume,  and  represents  the  equilibrium  pressure 
of  an  ideal  gas  undergoing  adiabatic  expansion  and  recompression. 
In  this  report,  various  "models"  of  the  explosion  bubble  interior 
are  examined  to  determine  the  extent  to  which  various  effects 
alter  the  overall  behavior,  with  the  objective  of  evaluating  the 
error  committed  by  adopting  one  or  another  of  these  models  for 
use  in  MACYL6  calculations.  In  all  cases,  the  water  outside  the 
bubble  is  considered  incompressible  as  ir.  the  MACYL6  cods,  and 
the  bubble  is  assumed  to  ramain  spherical  throughout  its  motion. 

A  bubble-interior  "model"  which  adequately  describes  a  spherical 
explosion  bubble  should  describe  the  more  general  case  equally 
well . 

Some  of  tha  explosion  bubble  models  which  will  be 
discussed  in  the  subsequent  sections  have  been  investigated  in 
the  past  by  other  workers  in  the  field.  The  "gasless"  case  dis¬ 
cussed  in  section  V  has  been  studied  by  many  authors,  among  them 
Lamb  (1932),  Willis  (1941)  and  Cole  (1948).  The  ideal  gas  model 
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of  section  VI  has  been  treated  by  Friedman  (1347)  and  by  Snay 
and  Christian  (1952)  at  some  length.  Many  of  these  early  results 
were,  however,  limited  in  precision  since  the  solutions  neces¬ 
sarily  involve  extensive  numerical  integrations .  Therefore,  the 
results  presented  herein  were  all  computed  quite  accurately 
using  very  finely-resolved  numerical  integration  procedures  on 
a  CDC-6600  computer.  Furthermore,  the  results  for  all  cases 
are  presented  in  a  consistent  way  to  facilitate  comparison  among 
the  various  models  of  bubble  behavior. 

In  the  final  section,  the  results  of  these  calculations 
are  compared  both  with  experimental  measurements  and  with  a  few 
compressible  calculations  of  the  water  motion  carried  out  by 
other  investigators;  in  general,  the  incompressible-water 
approximation  is  seen  to  be  good.  Before  proceeding  to  the 
development  of  the  various  bubble  models,  however,  the  general 
phenomenology  of  underwater  explosions  will  be  qualitatively 
described. 
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II . 


FORMATION  OF  THE  EXPLOSION  BUBBLE 


The  detonation  of  a.  submerged  explosive  sets  into 
motion  a  complex  sequence  of  events  of  which  the  character 
depends  on  the  nature  of  the  explosive,  the  energy  released, 
the  depth  of  burst,  and  other  environmental  parameters.  Nuclear 
explosions  are  of  primary  interest  in  this  study,  but  mc>3t  of 
the  experimental  information  available  concerning  underwater 
explosion  effects  was  gathered  using  chemical  explosives.  There¬ 
fore,  in  this  discussion,  both  types  of  explosions  will  be 
described,  taking  note  of  the  differences  in  effects. 

First,  we  will  consider  an  uncased  spherical  charge  of 
conventional  high-explosive  (ouch  as  TNT)  initiated  at  the  center. 
As  the  detonation  front  expands  through  the  charge,  the  solid 
explosive  encompassed  undergoes  chemical  reaction  and  releases 
energy  to  further  drive  the  detonation  shock.  The  detonation 
wave  speed  is  typically  in  the  range  6000-7000  meters  per 
second  for  most  explosives.  Behind  the  detonation  front,  the 
"burned"  reaction  products  will  include  such  materials  as  CO, 

C02,  H20,  NO,  CH^  and  «2  as  gases,  and  C,  Pb,  and  A1203  as 
solids.  Ultimately,  of  course,  the  detonation  wave  will  reach 
the  charge  surface  and  proceed  into  the  water  as  a  strong  hydro- 
dynamic  shock  wave  at  which  time  the  chemical  reactions  will  be 
complete.  For  a  given  explosive  type  and  packing  density,  the 
total  energy  released  by  the  explosion  will  be  proportional  to 
the  charge  mass;  that  ist 

Y  -  y  s  R*  pB  0  (II-l) 


where 

Y  ■  yield;  total  explosion  energy 
Rq  -  charge  radius 
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P_  *  explosive  packing  density 

h> 


Q  «  chemical  energy  released  per  unit  explosive 
mass . 

Furthermore,  the  "reaction  time"  (time  required  for  the  detona¬ 
tion  process)  will  depend  on  the  charge  size  and  th».  detonation 
wave  speed; 


rJ>_  . 

'det 


3Y 


^E^det 


1/3 


(II-2) 


For  TNT,  for  example,  a  typical  packing  density  is  about  1500 

kg/m^  (that  is,  a  specific  gravity  of  1.5)  and  Q  is  approximately 

4.2  x  io6  joules/kg.  Therefore,  a  spherical  TNT  charge  with  an 

12 

energy  release  of  one  kiloton  (defined  as  10  calories  or  about 
12 

4.2  x  io  joules)  would  be  about  10.9  meters  in  diameter,  and  the 
total  reaction  time  would  be  slightly  less  than  one  milli«econd 
(see  Figure  1). 


Initially,  the  shock  propagated  into  the  water  contains 
about  half  the  explosion  energy;  the  remainder  resides  as  both 
kinetic  energy  associated  with  the  expansion  of  the  gaseous 
"bubble"  of  reaction  products,  and  as  internal  energy  (heat) 
within  the  explosion  products  themselves.  This  shock  travels 
away  rapidly,  declining  in  strength.  Near  the  original  charge 
position,  the  highly  non-linear  behavior  of  the  shock  causes 
energy  dissipation  as  heat  to  the  water,  once  the  shock  is 
about  10  or  15  charge  radii  from  the  origin,  however,  this 
dissipation  process  is  largely  complete;  the  shock  thereafter 
continues  to  travel  away  carrying  with  it  about  one-fourth  of 
the  original  explosion  energy  (see  Figure  2) .  The  time  interval 
from  the  moment  of  explosion  to  the  end  of  this  "dissipation" 
phase  (defined  for  our  purposes  as  the  moment  when  the  shock 
has  propagated  10  charge  radii)  is  still  fairly  short,  but  is 
about  25  times  greater  than  the  "reaction  time". 
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SHOCKWAVE  ENERGY  /  TOTAL  YIELD 


TIME  (millisecond* ) 


FIGURE 


2:  LARLY  SHOCKWAVE  ENERGY  DISSIPATION  FOR  A  QNE-KILQTON 

TNT  CHARGE  DETONATED  UNDERWATER 


The  difference  between  the  effects  of  a  chemical 
underwater  explosion,  such  as  described  above,  and  those  of  a 
nuclear  or  thermonuclear  explosion  arise  from  both  the  (usually) 
higher  yield  in  the  nuclear  case  ar.d  the  far  higher  initial 
explosion  energy  density.  The  nuclear  reactions  may  be  regarded 
as  essentially  instantaneous;  the  materials  which  comprise  the 
nuclear  device  (largely  metals,  along  with  the  reaction  pro¬ 
ducts  of  the  high-explosive  stage,  any  unbumod  nuclear  fuel, 
and  the  radioactive  nuclear  reaction  products)  are  raised  to 
extremely  high  temperatures  and  ionized.  The  earliest  energy 
transfer  mechanism  to  the  surrounding  water  mass  is  therefore 
radiation;  since  water  is  opaque  to  bomb- temperature  photons 
(the  mean  free  path  is  only  a  centimeter  or  so) ,  this  radiation 
transfer  process  may  be  adequately  described  as  "radiation 
diffusion".  The  resulting  high  pressurs  region  then  forms  an 
extremely  strong  hydrodynamic  shock  wave  which  expands,  encom¬ 
passing  more  and  more  water  and  raising  Its  entropy.  Near  the 
burst  point,  the  internal  energy  increase  per  unit  mass  imparted 
to  the  water  is  sufficient  that,  upon  expansion,  the  water  will 
vaporize.  This  internal  energy  jump  declines  as  the  shook 
strength  decreases  due  to  its  increase  in  surface  area,  however, 
and  therefore  the  shock  front  energy  density  eventually  becomes 
low  enough  that  no  further  water  will  be  vaporised.  It  turns 
out  that  this  separation  between  the  shook  front  and  the  bubble 
front  occurs  at  a  radius  roughly  equal  to  that  of  a  Spherical 
TNT  charge  of  the  same  energy  release  es  the  nuclear  explosion. 
Thus,  at  this  stage,  we  havs  a  physical  situation  not  entirely 
unlike  that  of  the  high  explosive  burst  at  the  moment  when  the 
detonation  wave  reeohes  the  charge  surface.  The  shock  wave  will 
thereafter  continue  to  propegete  outward,  diseipating  energy  ae 
haat  at  olose-in  distances,  and  rapidly  becoming  weaker,  similar 
to  the  high-exploelve  ease.  A  nuclear  explosion  leaves  about 
401  of  ita  anergy  behind  as  "initial  bubble  energy";  moat  of  the 
remainder  is  shockwave  energy  (of  which  a  portion  is  dissipated 
near  the  burst  point)  and  a  relatively  smell  fraction  remains 
as  energy  to  be  released  later  by  the  decay  of  the  radioactive 
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fission  products.  This  may  be  compared  with  the  47-53%  energy 
partition  between  the  "bubble"  and  the  primary  shockwave  for 
TNT  (Cole,  1949) ;  for  most  conventional  high-explosivea ,  the 
bubble  energy  fraction  is  between  40%  and  60%.  The  difference 
between  the  high-explosive  and  nuclear  cases  is,  of  course,  that 
the  interior  "atmosphere"  of  the  bubble  produced  by  a  conven¬ 
tional  underwater  explosion  is  initially  roughly  homogeneous 
and  is  composed  of  the  gaseous  reaction  products  of  the  charge, 
whereas  in  the  nuclear  case  the  bubble  atmosphere  consists  of 
steam  and  is  non-homoyeneous ,  being  much  hotter  at  the  center 
than  at  the  periphery. 


III. 


DESCRIPTION  OF  BUBBLE  PULSATION 


By  the  end  of  the  dissipation  phase,  the  shockwave 
has  propagated  well  away  from  the  "initial  bubble"  which  con¬ 
sists  in  the  conventional  high  explosive  case  of  the  charge's 
reaction  products  and  in  the  nuclear  case  of  the  "potential 
steam”  discussed  above.  After  this  point,  there  is  a  marked 
qualitative  similarity  between  the  two  cases.  If  the  explo¬ 
sion  is  sufficiently  shallow,  of  course,  the  bubble's  expansion 
will  rupture  the  surface,  throwing  up  a  hollow  vertical  column 
of  water,  and  violently  expelling  the  bubble  contents  into 
the  atmosphere.  A  typical  example  of  this  sort  of  explosion 
is  the  familiar  CROSSROADS-BAKER  nuclear  test  of  1946.  For 
deeper  explosions,  however,  the  bubble  will  not  vent  into  the 
air,  but  will  continue  to  expand,  at  first  rapidly,  and  then 
more  and  more  slowly.  The  internal  pressure  and  temperature 
will,  of  course,  drop  during  this  expansion,  and  at  some  point 
the  average  internal  pressure  will  become  equal  to  the  ambient 
hydrostatic  pressure  at  the  burst  depth.  The  momentum  of  the 
watc*  rushing  away  from  the  point  of  burst  will,  however,  carry 
the  expansion  even  further.  Eventually,  the  expansion  will  be 
brought  to  a  halt  by  the  hydrostatic  pressure,  and  the  bubble 
will  begin  to  collapse.  The  interface  will  move  inward  with 
ever-increasing  speed  and  recompression  of  the  bubble  atmosphere 
will  occur  until  the  motion  is  cnce  again  brought  to  a  halt  by 
the  high  internal  pressure;  thereafter,  the  bubble  will  re¬ 
expand  and  th*  '•’hole  process  will  be  repeated.  The  reversal  of 
the  motion  at  the  bubble  "minimum 9  is  so  abrupt  as  to  appear 
discontinuous  on  a  time  scale  appropriate  for  the  expansion- 
contraction  cycle  as  a  whole.  On  the  other  hand,  the  motion 
near  the  "maximum"  is,  relatively  speaking,  extremely  smooth 
and  leisurely;  the  bubble  radius  is  more  than  half  the  value 
at  the  bubble  maximum  for  over  90%  of  the  period  of  oscillation. 

During  most  of  the  expansion-contraction  cycle,  there 
is  little  opportunity  for  energy  exchange  between  the  bubble 
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interior  and  the  surrounding  water.  Over  the  time  periods  in 
question,  heat  conduction  across  the  bubble  boundary  will  be 
entirely  negligible,  and  hence  the  bubble  atmosphere's  expansion 
and  recompression  may  be  taken  as  adiabatic.  Furthermore,  a 
one-kiloton  nuclear  explosion  at  a  burst  depth  of  100  meters, 
for  example,  will  generate  a  steam  bubble  whose  radius,  at 
maximum  expansion,  is  about  65  meters  and  whose  period  of 
oscillation  is  about  3.8  seconds.  Thus,  it  can  be  seen  that 
the  bubble-pulsation  effect  iB  a  relatively  slow-motion 
phenomenon  and  is  confined  to  a  fairly  small  region  of 
space,  compared  to  the  more  familiar  effects  of  the  shock  wave. 

An  overall  "characteristic  velocity"  for  the  oscillation  may 
be  taken  as  simply  the  maximum  radius  divided  by  the  period  of 
oscillation,  or  about  17  meters/second  in  the  above  typical  case. 
This  is  two  orders  of  magnitude  smaller  than  the  speed  of  sound 
in  water,  and  consequently  the  water  motion  may  be  adequately 
treated  as  incompressible  over  most  of  the  bubble  cycle. 

Near  '.he  bubble  minimum,  on  the  other  hand,  the 
assumptions  of  water  incompressibility  and  adiabatic  gas  behavior 
begin  to  break  down.  First,  very  close  to  the  minimum,  the  water 
velocities  adjacent  to  the  bubble  interface  become  very  high, 
and  u  weak  pressure  pUj.*u  iv  rauxfirvo  away  from  tne  bubble , 
carrying  with  it  a  few  percent  of  che  bubble's  pulsation  energy. 
This  "bubble  pulse"  is  somewhat  broader  but  is  much  lower  in 
amplitude  than  the  primary  shockwave.  Second  (and  much  more 
important) ,  at  the  minimum,  Taylor  instability  occurs  at  the 
water-bubble  interface.  This  instability  is  that  of  capillary 
waves  on  the  bubble  surface i  thus,  the  size  of  the  perturbations 
is  quite  small  compared  to  the  bubble  size  even  at  its  minim'im, 
and  furthermore,  the  instability  is  a  fully  three-dimensional 
phenomenon.  The  resulting  interface  breakdown  causes  the  forma¬ 
tion  of  a  spray  of  water  droplets  which  penetrate  the  bubble, 
cooling  its  interior.  If  the  explosion  is  nuclear,  this  cooling 
will  cause  condensation  of  a  portion  of  the  steam  atmosphere. 


thereby  causing  an  additional  energy  loss.  Finally,  the  in¬ 
stability  of  the  interface  initiates  the  development  of  intense 
turbulence  adjacent  to  the  bubble,  and  the  energy  of  this  turbu¬ 
lence  is,  of  course,  derived  from  bubble  energy.  Thus,  the  second 
cycle  of  oscillation  is  weaker  than  the  first,  and  each  subsequent 
oscillation  is  damped  even  further.  Experimental  measurements 
have  shown  that  if  the  bubble  does  not  migrate  upward  appreciably 
due  to  gravity,  the  bubble  energy  available  for  the  second  cycle 
of  oscillation  is  about  40%  of  that  of  the  first  in  the  high- 
explosive  case;  the  second  cycle  energy  of  a  steam  bubble, 
with  its  condensible  atmosphere,  is  only  about  8%  of  the  first 
cycle  energy  under  the  same  circumstances.  In  the  high-explosive 
case,  the  bubble  may  experience  as  many  as  eight  or  ten  oscilla¬ 
tions  before  becoming  relatively  inert;  in  the  nuclear  case, 
however,  the  steam  bubble  will  generally  condense  awiy  completely 
at  tha  end  of  the  third  cycle.  Thif  behavior  is  illustrated  in  ; 
Figure  3.  -  -  - 


FIGURE  3 i  THE  EFFECTS  OF  BUBBLE  ENERGY  LOSS 


NON-MIGRATING  CASE 


IV. 


TWO-DIMENSIONAL  EFFECTS 


So  far,  it  has  been  tacitly  assumed  that  the  bubble 
remains  spherical  in  form  and  that  its  center  does  not  move 
significantly  during  the  motion.  For  most  cases  of  interest, 
however,  these  assumptions  are  not  entirely  valid.  As  has  been 
mentioned,  if  the  explosion  is  shallow  enough,  the  bubble  may 
rupture  the  surface  early  in  its  expansion  and  expel  its  con¬ 
tents  into  the  atmosphere.  Furthermore,  even  if  the  explosion 
is  well  underwater,  if  it  occurs  sufficiently  close  to  (or, 
particularly,  in  contact  with)  the  sea  bed,  the  motion  would 
be  expected  to  be  influenced  by  this  solid  boundary.  Therefore, 
in  either  of  these  cases,  the  motion  is  not  one-dimensional 
(that  is,  the  state  of  the  system  is  not  definable  at  a  parti¬ 
cular  txme  in  terms  of  the  distance  from  the  burst  point  alone) 
but  is  axisymmetric  (requiring  two  space  coordinates,  r  in  the 
radial  direction  from  the  vertical  axis  of  symmetry,  and  z, 
the  altitude) . 

Even  if  the  explosion  is  deep  enough  that  the  surface 
is  not  ruptured,  and  far  enough  from  the  sea-floor  that  the 
presence  of  this  boundary  does  not  significantly  influence  the 
flow,  the  bubble  motion,  in  general,  will  still  not  be  one¬ 
dimensional.  The  reason  for  this  is  the  effect  of  gravity;  that 
is,  the  buoyancy  of  the  bubble.  As  a  consequence  of  this  buoy¬ 
ancy,  the  bubble  will  tend  to  float  upward  toward  the  wat?r 
surface  as  it  pulsates.  The  instantaneous  buoyant  force  is 
given  by  Archimedes'  principle: 

Fb  “  (pw  ~  pg)gV  UV-D 


where 


Pw  ■  density  of  water  (constant) 
pg  ■  density  of  bubble  interior 
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g  =  acceleration  of  gravity  (constant) 

V  =  bubble  volume. 

The  upward  momentum  generated  by  buoyancy  up  to  a  time  t*  is 
therefore  simply: 

ft* 

I  =  g  |  <pw  -  p  )V  dt  ( IV-2 ) 

O'  ’ 

During  most  of  the  bubble  cycle,  when  V  is  large,  p^  (the  gas 
density)  will  be  orders  of  magnitude  smaller  than  pw>  There¬ 
fore,  for  our  purposes,  we  may  say: 

ft* 

I  =  p  g  |  V  dt  (IV- 3) 

This  momentum  is  accumulated  with  time  as  the  motion  proceeds. 
Near  bubble  maxima,  when  V  is  large,  momentum  is  accumulated 
most  rapidly.  The  "hydrodynamic  mass"  or  "mass  of  moving 
water"  is,  however,  also  very  near  maxima.  At  the  Inst  minimum, 
on  the  other  hand,  the  "hydrodynamic  mass"  is  small,  but 
possesses  all  the  upward  momentum  accumulated  during  the  first 
cycle,  that  is, 

I  -  Pw9  I"  V  dt  =  |  7rpwg  [T  R*dt  ( IV-4 ) 

o  o' 

where  T  is  the  first  osc  llation  period  and  R  is  the  bubble 
radius.  Therefore,  the  overall  upward  velocity  is  greatest  at 
bubble  minima. 

The  volume-time  integral  for  the  first  bubble  cycle 
in  equation  (IV-4)  may  be  taken  as  simply  proportional  to  the 
volume  at  the  maximum  times  the  oscillation  period,  for  a  pre¬ 
scribed  radius-time  curve: 

fT 

R’dt  -  B  Rj*^  T  (IV-5) 

o' 


22 


whore  1^  is  the  first  maximum  bubble  radius,  and  the  coeffi¬ 
cient  of  proportionality  B  depends  upon  the  exact  shape  of  the 
radius-time  relation  for  the  bubble  oscillation.  If  a  "charac¬ 
teristic  momentum"  for  the  bubble  oscillation  is  defined  as: 


4 

7 


rr 


(IV-6) 


the  upward  "buoyant"  momentum  may  be  normalized  and  presented 
in  dimensionless  form: 

I*  -  f2^-  (IV- 7) 

max 

which  will  be  recognized  as  a  reciprocal  Froude  number  (increasing 
with  increasing  relative  buoyancy  effect).  That  is,  two  explo¬ 
sions  which  produce  bubbles  such  that  that  quantity  I*  is  the 
same  will  experience  the  same  relative  upward  migration  due  to 
gravity  during  the  first  cycle.  Clearly,  therefore,  the  quantity 
B  (which  is  a  measure  of  the  shape  of  the  radius-time  curve)  is 
an  important  parameter  in  the  study  of  underwater  explosion  bubble 
pulsation. 


It  has  long  been  recognized  experimentally  that  if 
migration  is  strong  (that  is,  if  I*  is  large) ,  the  bubble, 
while  initially  spherical,  will  become  more  and  more  non- 
spherical  as  the  motion  proceeds.  After  the  first  maximum,  the 
bubble  bottom  will  tend  to  collapse  back  toward  the  explosion 
point  more  quickly  than  does  the  top,  thereby  generating  an 
upward  central  jet  of  water  which  collides  with  the  bubble  top 
just  prior  to  the  moment  of  maximum  recontpression.  Therefore, 
at  the  minimum,  the  bubble  is  toroidal  rather  than  spherical  in 
form.  If  the  migration  strength  is  sufficiently  great,  the 
bubble  may  remain  toroidal  thereafter;  otherwise,  the  central 
jet  will  dissipate  upon  re-expansion  but  may  re-form  at  sub¬ 
sequent  minima.  For  very  strongly  migrating  bubbles,  the  net 
upward  displacement  by  the  end  of  the  first  bubble  cycle  may 
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be  comparable  to  the  first  maximum  bubble  radius. 

The  MACYL6  hydrodynamic  code  (Pritchett,  1970a) ,  although 
fairly  general  in  application,  was  specifically  designed  to  deter¬ 
mine  the  bubble  motion  and  associated  phenomena  following  deep 
underwater  explosions.  The  MACYL6  code  is  a  computer  program  fo- 
solving  viscous  turbulent  time-dependent  incompressible  axisym- 
metric  fluid  flows  involving  free  surfaces.  The  fundamental 
ensemble-averaged  Navier-Stokes  equations  which  govern  all  such 
flows  are  numerically  integrated  using  finite-difference  methods 
on  an  Eulerian  mesh  of  computational  grid  points.  Free  surfaces 
are  treated  using  the  MAC  ( "Marker -and-Cell")  technique  first 
developed  by  Welch,  Harlow,  Shannon  and  Daly  (1966)  at  Los  Alamos. 
In  this  procedure,  the  fluid  is  "tagged"  with  a  large  number  of 
massless  "marker  particles"  which  are  moved  with  the  flow  through 
the  Eulerian  mesh  at  each  time  step  using  velocities  interpolated 
from  nearby  principal  grid  points.  These  marker  particles  there¬ 
by  delineate  the  positions  of  the  free  surfaces.  A  heuristic 
model  to  determine  the  effects  of  turbulence  was  developed 
separately  (Gawain  and  Pritchett,  1970)  and  is  an  integral  part 
of  the  overall  procedure.  For  application  to  the  underwater 
nuclear  explosion  problem,  the  fluid  is  taken  as  initially  at 
rest  and  the  water  surface  is  horizontal,  but  one  or  two  "empty" 
cells  at  the  explosion  point  contain  a  high  internal  pressure. 

The  cof,e  then  calculates  the  subsequent  water  flow  in  a  stepwise 
fashion  using  the  f inite-dif ferenc4?  forms  of  the  governing  equa¬ 
tions  and,  as  boundary  conditions,  the  fact  that  the  air  pressure 
is  constant  and  that  the  bubble  internal  pressure  is  a  prescribed 
function  of  bubble  volume.  For  all  such  calculations  published 
to  date,  the  bubble  pressure  -  volume  relation  for  a  particular 
cycle  of  oscillation  was  assumed  to  be  that  of  the  adiabatic 
rarefraction  and  recompression  of  an  ideal  gas  with  y  (the  ratio 
of  specific  heats)  equal  co  4/3,  that  is, 

PVV3 


»  constant 


(IV-8) 


As  was  discussed  in  section  III,  although  the  bubble's 
expansion  and  recompression  may  be  taken  as  adiabatic  during 
most  of  the  oscillation  cycle,  energy  is  lost  near  the  bubble 
minimum  due  to  two  fundamental  causes.  First,  near  the  minimum, 
a  weak  pressure  wave  (often  called  the  "bubble  pulse"  to  dis¬ 
tinguish  it  from  the  primary  shockwave)  is  emitted  from  the 
bubble  and  carries  away  a  small  fraction  of  the  bubble  energy. 

This  occi  rs  due  to  the  compressibility  of  water;  near  the 
minimum,  the  bubble  interface  velocity  may  instantaneously  become 
comparable  with  the  speed  of  sound  in  water,  even  though  over  the 
remainder  of  the  pulsation  velocities  are  far  below  sound  speed. 
Second,  the  reversal  of  the  motion  at  the  minimum  causes  Taylor 
instability  of  the  bubble-water  interface,  which  generates  intense 
turbulence  in  the  water  adjacent  to  the  bubble  and  causes  the 
breakdown  of  the  interface  into  spray,  which  in  turn  cools  the 
bubble  interior  and  (for  nuclear  explosions)  causes  condensation 
of  a  fraction  of  the  steam  interior.  This  instability  is,  how¬ 
ever,  a  fully  three-dimensional  phenomenon  and  the  scale  of  size 
of  the  perturbations  (i.e.,  the  spray  droplets)  is  far  smaller 
than  the  bubble  itself.  It  therefore  seems  clear  that  an  incom¬ 
pressible  spherically  symmetric  or  axiaymmetric  treatment  of  the 
water  motion  cannot  possibly  predict  the  energy  loss  at  the  bubble 
minimum.  Fortunately,  there  is  no  particular  requirement  for  such 
predictions;  the  fractions  of  the  bubble  energy  lost  at  each 
minimum  have  been  determined  experimentally  and  have  baen  presented 
by  Phillips  and  Snay  (1968)  for  steam  bubbles  and  by  Snay  (1962) 
for  tne  conventional  high-explosive  case,  It  turns  out  that  the 
fraction  of  the  bubble  energy  lost  at  the  bubble  minimum  declines, 
in  general,  with  increasing  migration  strength.  In  the  MACYL6 
code,  the  constant  in  equation  (IV-8)  is  adjusted  downward  at  each 
bubble  minimum  in  such  a  way  that  the  fractions  of  tha  bubble 
energy  lost  at  each  minimum  correspond  to  Snay's  empirical  results; 
thus,  the  decay  of  the  bubble  oscillation  with  time  is  properly 
reproduced  in  the  calculations. 
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An  earlier  version  of  this  code  has  been  used  with 
considerable  success  in  the  calculation  of  the  bubble  motion 
after  the  WIGWAM  deep  underwater  nuclear  explosion  test  of 
1955  (Pritchett  and  Pestaner,  1969)  .  The  MAC  If  L  6  code  itself  has 
been  used  more  recently  tc  compute  the  bubble  motion  and  subsequent 
flow  after  extremely  deep  underwater  nuclear  bursts  out  to  several 
minutes  after  detonation  (Pritchett ,  1971)  and  is  now  being  used 
to  study  the  formation  of  surface  waves  by  relatively  shallow 
Nuclear  explosions.  As  was  seen  above,  however,  the  MACYL6  code 
actually  computes  only  the  water  flow  outside  the  explosion 
bubble;  the  gas  dynamics  of  the  bubble  interior  are  represented 
only  in  an  approximate  way  (such  as  by  equation  IV-8) .  Although 
this  procedure  has  been  quite  successful  at  accurately  computing 
bubble  motion  (Pritchett .  1970b),  it  is  certainly  worthwhile  to 
investigate  the  adequacy  of  assumptions  such  as  (IV-8) ,  so  as  to 
ascertain  whether  or  not  an  improved  "bubble  interior  model"  can 
be  formulated  which  will  significantly  improve  the  results.  For 
this  purpose,  however,  it  is  not  necessary  to  make  fully  two- 
dimensional  (that  is,  axisymmetric)  calculations.  If  the  one- 
dimensional  case  can  be  adequately  treated  using  a  particular 
"model"  for  the  bubble  interior,  the  two-dimensional  case  in  which 
the  bubble  may  become  non-spherical  will  also  be  adequately  des¬ 
cribed  by  the  same  sort  of  model. 

Therefore,  we  will  consider  the  case  of  an  explosion  in 
an  unbounded  mass  of  water  in  which  the  ambient  hydrostatic  pressure 
is  everywhere  the  same.  The  effects  of  viscosity  (which,  in 
reality,  are  exceedingly  small)  will  be  neglected,  and  furthermore, 
the  water  outside  the  bubble  will  be  considered  incompressible. 

As  was  inferred  in  the  previous  discussion,  thin  assumption  is 
valid  except  for  very  brief  time  intervals  near  the  bubble  minima. 

In  any  case,  the  objective  of  this  study  is  to  determine  optimum 
procedures  for  representing  the  effects  of  the  bubble  interior 
in  the  MACYL6  code,  and  MACYL6  treats  the  water  in  the  incompres¬ 
sible  approximation. 
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V. 


SOLUTIONS  NEGLECTING  THE  BUBBLE  ATMOSPHERE 


The  first  (and  most  primitive)  model  for  the  bubble 
interior  which  will  be  considered  iB  that  there  is  no  bubble 
atmosphere  at  all;  that  is,  that  the  motion  begins  with  water 
rushing  away  irom  an  infinitesimal  point  forming  an  evacuated 
cavity.  To  derive  the  governing  equations  for  the  resulting 
motion,  we  first  impose  Euler's  momentum  conservation  equation 
for  an  incompressible  fluid  in  spherically  symmetric  motion: 


3u  1_  3 

3t  *  Tr 


(r2u2 ) 
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3P 

TI 


where 

u  •»  radial  velocity 
r  *  distance  from  origin 
P  “  pressure 

Pw  "  fluid  density  (constant) ,  and 
t  -  time. 


(v-l) 


Furthermore,  the  continuity  condition  is: 

^7  37  (r!u>  '  0  (V-J> 

If  at  a  particular  time  the  radius  of  the  babble  is 

■ 

denoted  by  R  and  its  instantaneous  velocity  of  expansion  is  R, 
the  velocity  elsewhere  in  the  water  is  just: 

u  *  R(R/r)2  for  r  >  R  (V-3) 


as  can  be  seen  from  the  above  continuity  condition.  Inserting 
this  result  into  the  momentum  equation,  we  immediately  obtain: 
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where 


p  .=  ,  *1* 

at2 

as  the  governing  equation. 


To  eliminate  pressure,  we  try  a  solution  of  the  form: 

(V-5) 


P_ 

P, 


al  a2  a3  a4 

°0  +  —  +  ~+  —  4-1 

0  r  rJ  r*  r1* 


so  that: 


1  3P 
Pw  3r 


«X  ^2  3a3  4a4 

ri  j-l  j.4  rS 


(V-6) 


Substituting  the  above  into  the  governing  equation  (V-4) ,  we 
obtain: 


a1  -  R  R2  +  2  R2R 


a2  "  a3  "  0 


«4  -  '  j  RJR“ 


(V-7) 


and  therefore : 


£_  -  ,  +  U!L  + 

Pw  o  r 


*D* 


2r*  5  -  El* 


2r* 


(V-8) 


Now,  at  r  -  »,  the  pressure  is  simply  the  hydrostatic 


pressure,  PH: 


and  thus: 


P  -  PH  at  r 


(V-S>) 
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I  i 

l: 

I 

I 

I 


(V-10) 


Therefore/  the  pressure  field  is  given  by: 


P 


PH  +  PW 


+  2R2 


H 

r 


£i£l) 

2r“ 


(V-ll) 


If  we  define  the  instantaneous  pressure  at  the  bubble-water 
interface  as  PB,  we  may  evaluate  the  above  equation  at  r  -  R 
to  obtain  an  equation  of  motion  for  the  bubble: 


R  = 


(V— 12) 


Now,  if  the  bubble  atmosphere  is  ignored,  the  bubble 
pressure  PQ  is  always  zero;  thus,  finally,  we  obtain: 


R 


3  R2 
7  R“ 


(V-13) 


It  is  also  feasible  to  derive  the  governing  equation 
from  the  energy  principle.  The  total  energy  in  the  system  re¬ 
mains  constant,  and  is  composed  of  two  terms;  the  kinetic  energy 
of  the  moving  water,  and  the  total  work  done  against  external 
forces  (in  this  case,  the  hydrostatic  pressure); 


eq  =  er  +  H  (V-14) 

The  total  kinetic  energy  may  be  evaluated  by 

Ek  -  \  1  pw  u2<iv  *  2lrPw  j  r2u2dr  fv~15) 

Using  the  continuity  condition  (equation  V-3  above) ,  this  becomes: 

ek  -  2"»„  *'*•  J"  ff 

-  2TrpwR2R*  (V-16) 
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The  work  done  against  hydrostatic  pressure  to  expand  the  bubble 
from  radius  RQ  to  radius  R  is  simply: 

fR 

W  =  j  Pj^  dR 

where  A  is  the  instantaneous  bubble  area,  or: 

/R 

W  =  4ttPh  R2dR 

«  |  trPH  (V-17) 

In  the  present  case,  the  bubble  starts  from  zero  radius,  and 
therefore  the  total  work  done  is  just: 

W  “  7  (V-18) 

Therefore,  the  energy  equation  becomes: 

Eq  *  2irpw  R2RS  +  j  uPjjR5  (V-19) 

If  this  equation  is  differentiated  with  respect  to  time,  the 
result  is: 


0  *»  4i rpw  R  R  R'  +  6ttpw 


R’R2  +  4ttPhR  R2 


which  may  be  solved  algebraically  for  R: 


f  p 

H 

*7 


.  3  ft2 

+  7  IT 


which  is  identical  to  the  governing  equation  derived  from  the 
momentum  principle  (equation  V-13). 


A  special  case  of  the  "gaslesa"  bubble  solution  is  that 
which  is  obtained  if  the  hydrostatic  pressure  is  ignored 7  that  is, 
in  terms  of  energy  equation,  the  kinetic  energy  remains  constant: 


zStJL  v  -  2$  ' 
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<V-20) 


Eo  “  2TrP„  r2r’ 


for  which  the  solution  may  be  obtained  analytically: 


R 


1/5 

X 


t2/5 


(V-21) 


in  this  case,  the  bubble  radius  increases  monotonically  with 
time,  and  no  bubble  pulsation  occurs  due  to  the  absence  of  hydro¬ 
static  pressure.  This  result  is  approximately  valid  at  early 
tiroes,  but  becomes  worse  and  worse  as  the  bubble  grows  and  the 
work  done  against  hydrostatic  pressure  becomes  an  important  term 
in  the  bubble  energy  balance. 


This  case  is  interesting,  however,  in  that  it  permits  a 
rough  estimate  of  the  adequacy  of  the  incompressible  assumption. 
The  error  may  be  estimated  as  being  of  the  order  of  the  square  of 
the  Mach  number  based  upon  the  interface  velocity,  which  may  be 
shown  from  (v  -21)  to  be: 


M' 


1 

57 


PwR  C 


(V-22) 


where  C  is  the  speed  of  sound  in  water.  This  quantity  is  plotted 
in  Figure  4  for  the  one-kiloton  case  discussed  earlier;  as  can 
be  seen,  at  a  radius  corresponding  to  the  TNT-equivalent  charge 
radius,  the  local  Mach  number  is  significantly  less  than  unity, 
similarly,  the  square  of  the  local  Mach  number  as  a  function  of 
time  is : 


M1 


4 

25C* 


.-6/5 


<V-23> 


which  is  plotted  for  a  one-kiloton  explosion  in  Figure  5.  The 
incompressible  approximation  is  clearly  wr -ranted  after  about  5 
or  10  milliseconds;  the  overall  oscillation  period  of  the  explo¬ 
sion  bubble  is  several  seconds.  Therefore,  these  results  tend  to 
support  the  conclusion  drawn  earlier  than  the  water  flow  around 
the  bubble  may  be  adequately  treated  as  incompressible  motion. 
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If  the  hydrostatic  pressure  is  non-zero,  of  course,  the 
bubble  will  not  expand  indefinitely,  but  will  oscillate  in  size. 
If  a  "reduced  radius”  and  "reduced  time"  are  defined  as  follows: 


P„  1/3 
A  «  R  (*— ) 

Eo 


(V-24) 


T  -  t 


(V-25) 


o  '  E 
w  o 


and  these  variables  substituted  into  equation  V-19,  the  energy 
equation  may  be  written  in  dimensionless  form: 


2n  A  2  A  *+  |  it  A1  -  1 


(V-25) 


where 


i  -  dA 
x  a? 

As  can  be  seen,  the  maximum  radius  of  the  bubble  will  occur  when 
the  bubble  interface  velocity  is  zero: 


A  -  (f-)  -  0.620 


(V-26) 


or,  in  dimensional  form. 


E  1/3 

■W  '  <*;> 


(V-27) 


where 


■  maximum  dimensionless  radius  A 
-  maximum  bubble  radius 

»  bubble  energy  (about  half  the  total  explosion 
energy) 

«  hydrostatic  pressure 


The  energy  equation  (V-  5)  may  be  integrated  numerically: 
the  resulting  dimensionless  radius-time  relation  is  illustrated  in 
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Figure  6.  The  relative  contributions  of  the  kinotic  and  internal 
energy  terms  to  the  total  energy  as  functions  of  time  are  shown 
in  Figure  7.  As  is  to  be  expected  from  the  form  of  equation  V-25, 
the  solution  is  symmetric  around  the  time  of  the  bubble  maximum 
and  there  exist  periodic  singular  discontinuities  corresponding  in 
time  to  the  bubble  minima.  The  period  of  the  oscillation  in  dimen¬ 
sionless  form  is: 


6t  -  1.135 


(V-28) 


or,  in  dimensional  form. 


T  -  1.135  o. 


1/2  1/3  5/6 


(V-29) 


where  6t  is  the  t  interval  between  successive  zeroes  of  A,  and  T 
is  the  corresponding  dimensional  period  of  oscillation. 

Another  quantity  of  interest  which  has  been  discussed 
previously  is  the  "characteristic  velocity"  associated  with  the 
bubble  oscillation.  This  may  be  taken  as: 


°c  •  w* 

which,  using  V-27  and  V-29,  becomes: 


(V-30) 


*  P»  i/2 

u  =  V~'  (— ) 
c  St  vp.  / 


(V-3H 


that  the  corresponding  dimensionless  parameter  is: 


y  «*  A/St  ■  0.547 


(V-32) 


Therefore,  the  "characteristic  velocity”  is  independent 
of  the  explosion  yield,  and  depends  only  on  the  hydrostatic  pressure, 
which  in  turn  depends  on  the  depth  of  the  explosion  (the  water 
density  pw  is  assumed  constant).  Numerically,  this  characteristic 
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FIGURE  6  s  THE  GASLESS  MODEL  -  THE  DIMENSIONLESS  BUBBLE  RADIUS 
TIME  RELATION 


T 


T 


FIGURE  7:  THE  GASLESS  MODEL  -  THE  BUBBLE  ENERGY  BALANCE 


velocity  is  fairly  low,  aa  discussed  in  section  III,  and  is  far 

below  the  speed  of  sound.  Typical  values  for  various  explosion 

3  2 

depths,  assuming  that  •  1000  kg/m  ,  g  ■  9.0  m/sec  ,  and  that 

w  5 

the  air  pressure  over  the  water  is  one  atmosphere  (1.013  x  10 
2 

newtons/m  )  are  listed  in  the  following  table: 


Burst  Depth 

Characteristic  Velocity 

Square  of  Over: 

(Meters) 

(Maters/Sec) 

Mach  Number 

10 

7.7 

2.7  x  io~5 

30 

10.9 

5.3  x  io~5 

100 

18.0 

1.4  x  10“4 

300 

30.2 

4.1  x  10‘  4 

1030 

54. S 

1.3  x  IQ-3 

3000 

94.1 

ro 

i 

O 

•-* 

X 

as 

a 

co 

Thus  it  seems  clear  that  water  compressibility  effects  are  not 
important,  even  for  explosions  at  very  great  depths. 

One  other  parameter  required  for  the  proper  description 
of  the  bubble  oscillation  is  the  quantity  B,  which  is  needed  to 
determine  the  upward  momentum  generated  by  the  bubble  due  to 
buoyancy,  h  is  defined  as  (see  equation  IV-5) : 

/T 

R'(t)dt  X(T)*dT 

o )  oJ _ 

B  "  Cx  T  A'4t 

The  upward  momentum  generated  by  the  bubble  during  the  first 
cycle  of  oscillation  in  normalised  form  is  (see  equation  XV-7) t 

and  the  value  of  B  obtained  from  the  integration  of  the  energy 
equation  is: 

B  -  0.625. 
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The  very  nature  of  the  governing  equation  for  thie 
relatively  simple  case  assures  that  the  radius-time  relation 
which  results  is  a  universal  function.  In  the  more  complicated 
treatments  considered  in  the  next  few  sections,  the  radius-time 
relation  will  be  found  to  depend  on  various  parameters,  so  that 
the  constants  which  characterize  the  radius-time  function,  that  is, 


A 

«  0.620 

(7-33) 

6t 

-  1.135 

(V-34) 

U 

*  0.547,  and 

(V— 3  5 ) 

B 

*  0.625 

(V-36 ) 

will  become  functions  of  those  parameters.  In  general,  however, 
the  same  system  of  nomenclature  will  be  retained,  so  that  results 
may  be  readily  compared  with  the  present  case. 
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VI.  THE  IDEAL  GAS  MODEL 

In  the  case  just  discussed,  the  internal  pressure  of 
the  bubble  was  approximated  as  zero.  The  bubble  pressure  is, 
in  fact,  low  compared  to  hydrostatic  pressure  during  most  of 
the  bubble  cycle.  None  the  less,  it  seems  clear  that  the 
solution  could  be  improved  if  some  plausible  method  were  used 
to  simulate  the  effects  of  the  gas  inside  the  bubble. 

In  particular,  if  we  begin  with  a  sphere  of  compressed 
ideal  gas  with  volume  VQ  and  with  the  fluid  initially  at  rest, 
the  gas  intornal  pressure  may  be  defined  by: 

P_vY  -  p  vY  -  constant  (VI-1) 

B  O  O 

that,  is,  an  adiabatic  expansion.  Here,  is  the  pressure 
associated  with  bubble  volume  V,  and  PQ  is  the  Initial  pressure; 
Y  (the  ratio  of  specific  heats)  characterizes  the  gas.  Further¬ 
more,  it  will  be  assumed  that  the  gas  remains  homogeneously  dis¬ 
tributed  within  the  bubble  throughout  the  motion,  and  therefore: 

u  -  R  £  for  r  <  R  (VI-2) 

As  before,  of  course,  the  velocity  field  outside  the  bubble  is: 

u  -  ft  (£) 1  for  r  >  R  (VI-3) 

If  the  initial  density  of  the  gas  is  denoted  by  p£,  then  the 
ga  density  inside  the  bubble  at  any  time  is: 

R 

Pg  -  PE  (j2) *  for  r  <  R  <Vl-4) 

and,  of  course,  the  watei  density  is: 

Pw  ■  constant  for  r  ;>  K  (VJ-5) 
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The  energy  balance  now  consists  of  four  components: 


EKG  +  EKW  +  Ej  +  w  ■  constant  !VI-6) 

that  is,  the  kinetic  energy  of  the  expanding  gas,  the  water 
kinetic  energy,  the  internal  (heat)  energy  of  the  bubble  atmos¬ 
phere  and  the  work  done  against  hydrostatic  pressure.  Two  of 
these  have  already  beer,  worked  out  (see  equations  V-15  through 
V-17)  : 

Exw  *  *">«  E‘»‘ 
w  ■  J  *PHtR’-K£> 

The  gas  kinetic  energy  is  just: 

EKG  '  7  J  >,  U’av 

bubble 
interior 

fR  , 

-  2np  —  r“dr 
9  «■  0J 

“  y  'npg 

and,  using  equation  (VI-4) 

ekg  -  I 

The  internal  energy  of  the  gas  is  simply: 

Ej  -  (VI-10) 

inserting  (VI-1) ,  this  may  be  written: 


(VI-7) 

(VI-8) 
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(VI-11) 


E,  - 


4* 

5TFTT 


lY 
,Ty=TT 


P  R_ 
o  o 


Initially,  the  system  is  at  rest  and  no  work  has  yet 
been  done  against  external  forces.  Therefore,  the  total  energy 
at  t  ■  0  (the  constant  on  the  right  in  equation  Vl-6)  is  just 
the  total  internal  energy  initially  present: 

\  m  T^TT  "  ekg  +  ekw  +  Ei  +  w  <V1‘12> 


The  bubble  oscillation  energy  EQ  is  just  the  amount  by  which 
the  initial  internal  energy  exceeds  the  internal  energy  required 
to  maintain  the  bubble  against  the  ambient  hydrostatic  pressure: 


E_  - 


PUV 
H  O 

TFTT 


<WVo 

(y-11 


4ir 

TTy^TT 


<Vph>r 


I 

o 


(VI-13) 


Thus,  for  example,  if  the  bubble  oscillation  energy  BQ  were 
zero,  the  initial  bubble  pressure  would  equal  hydrostatic 
pressure  and  no  motion  would  occur.  The  energy  equation  may 
now  be  expressed  explicitly  be  assembling  equations  VI- 7,  -8, 
-9,  -11,  -12,  and  -13: 


Eo  -  2tpw  R’R»  ♦  f  wp,  R*R’  ♦ 


P  R*^ 
4  ir  poRo 


Tv^r 


+  7  *ph  (R'-  po>  -  TT^IT  Vo  ,VI-14’ 


which  may  be  compared  with  the  energy  equation  for  the  previous 
"gasleas"  case  (V-19) .  The  corresponding  momentum  equation  ie 
simply  (compare  to  V-13) : 
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(VI-15) 


»  .  1  PQ,  _  PB“PH  3  R* 

B(1  +  ^  V  “P^r  -  7  r* 


Making  use  of  (VI-1)  and  (VI-4) ,  this  becomes 


,  .  R  >Y 

i  ^  <Vir3 

\  Pn 

+  k  ^ '<r>  *! 


-  PH>  -  I  *a]  / 


(VI-16) 


The  terms  involving  p£  reflect  the  inertial  effects 
of  the  gas  inside  the  bubble.  For  example,  the  ratio  of  the 
gas  kinetic  energy  to  that  of  the  water  is: 


"KG 

"KW 


I  (!o. , 

5  Pw  (R  } 


(VI-17) 


For  nuclear  explosions,  of  course,  the  initial  bubble  density 
(P_)  would  be  taken  as  the  same  as  that  of  water  (P  ) .  Even 

E  W 

for  high  explosives,  {Pt/P„l  will  rarely  exceed  1.6  or  so. 
Furthermore,  over  most  of  the  bubble  oscillation,  R  is  vastly 
greater  than  Rq.  Therefore,  it  would  seem  warranted  at  this 
stage  to  neglect  gas  inertial  effects;  the  effect  of  this 
approximation  will  be  examined  later  in  this  section. 


If  all  terms  involving  tue  gas  density  are  neglected, 
the  energy  equation  becomes: 

R„  *  (7-1) 

_  .  h‘r'  +  2-  n  r ,-2—  fp  i 

o  w 


rI 

E_  »  2*P„  }  *  lySr  (po  <1T> 


-  PH> 


+  P„ (R*  -  R* ) ] 
H  o 


(VI-18) 


and,  for  the  momentum  equation; 


R  *  4  “  4 

w 


»Y 


R  '‘o  »JTJ  rH'  t 


3 


(VI-19) 
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If  the  governing  equation  is  non-dimenaionalised  in 
the  same  way  as  waa  done  in  the  previous  section  for  the 
evacuated  bubble  case  by  the  introduction  of  the  dimensionless 
bubble  radius  (X)  and  the  dimensionless  time  (t),  defined  by: 

E  1/3 
R  -  X  U~) 
rH 


t  '  '  ^  ft 7S 

the  energy  equation  becomes : 


2ttX  2  X  s  +  3  ir{X  8 II  +  (^) 

-  -JL.  x8>-  i 

Y~1  O 


;vi-2o) 


which  may  be  compared  to  equation  (V-25)  for  the  "gasiess'’ 
case.  The  bubble  energy  is  just: 


_  OHO_«_  OH«) 

E  m  - - — - -  m  Tf  in— ■■■■—  ■  R 


(VI-21) 


or,  in  dimensionless  form. 


X*  - 


(y-l) 

Fw 


(VI-22) 


which  may  be  used  to  eliminate  XQ  from  the  energy  equation.  To 
facilitate  presentation,  we  define: 


I  ffPHRo 


(VI-23) 


which,  for  the  ideal  gas  case  under  discussion,  is  simply: 


P  "P„ 
«  OF 

c  “  TFITl^ 


(VI-24) 


i'hs  physical  significance  of  £  is  that  it  is  the  racio  of  the 
buuble  energy  to  the  work  required  to  form  the  initial  cavity 
of  radius  Rc  against  hydrostatic  pressure  Pj,.  Thus,  £  may  be 
leyai  Ced  as  a  measure  of  the  "relative  intensify"  of  the 
explosion,  or,  feu  a  given  explosive  type,  as  inversely  pro¬ 
portional  to  tlu  hydrostatic  pressure,  which  in  turn  varies 
vith  the  burst  depth.  In  particular,  the  range  of  interest  of  £ 
car.  be  estimated  by  using  TNT  as  a  standard  explosive  and  taking 

note  of  the  fact  that  the  initial  bubble  energy  density  for  TNT 
9  3 

is  about  3  *  10  joules/m  .  Thereby,  the  following  values  of  £ 
may  be  correlated  (approximately)  with  explosion  depth  as 
follows  t 


1 

I  1 
I 
1 
1 
1 


£ 

Depth  (meters) 

L  i 

1 

100 

3000 

1 .  i 

»  j 

300 

1000 

i 

1000 

300 

* 

2000 

90 

‘  i 

i 

10000 

20 

1 

1 

Actually,  as  will  be  ceen  later  on, 

eve  higher  values  oi  £ 

1 

i 

are  of  some 

interest,  since  much  of 

the  experimental  information 

!  ■  ! 

concerning 

bubble  puisati^  v  :s  collected  in  laberatory-scaie 

• 

test  chambers  in  which  a  vac  was 

drawn  over  the  water  to 

J| 

further  reduce  hydroat.  :ic  assure 

(see  for  example,  Buntze.., 

1 9 1  4  j  5 nay , 

1964  ;  Pritchett,  1S66)  . 

Thus,  the  range  of  variation 

r 

in  £  could 

extend  from  about  100  to 

as  high  as  100,000. 

L 

By  use  or  \Jl-22)  and  (VI-24),  the  dimensionless  energy 
equation  (Vi-20)  may  be  written  as  follows: 


2*a*a»  +  J  4  -J - '  TT^TP 


-  1 


(VI-24) 


As  can  bf*  seen,  the  solution  is  no  longer  a  universal  function 
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as  was  true  in  the  "gasless"  case  discussed  in  section  V,  but  is 
a  family  of  functions  which  depend  on  the  values  of  the  charac¬ 
teristic  parameters  y  (which  characterizes  the  gas)  and  f  (which 
can  roughly  be  correlated  with  explosion  depth) . 

Equation  (VI-24)  was  numerically  integrated,  using  a 
large  digital  computer,  for  values  of  y  varying  between  1,1  and  1.5, 
and  for  c  ranging  from  100  to  100,000.  Qualitatively,  the  radius¬ 
time  relations  obtained  were  quite  similar  to  that  of  the  "gasleHS" 
case  discussed  in  section  V;  the  result  for  y  *  4/3  and  c  -  3000 
is  illustrated  in  Figure  8.  The  relative  contributions  of  the 
various  terms  in  the  energy  equation  as  functions  of  time  for  this 
caaa  are  shown  in  Figure  9.  The  principal  qualitative  distinction 
between  the  two  cases  is  that  the  bubble  radius  doeB  not  drop  to 
zero  at  minima.  Quantitatively,  however,  the  constants  which 
characterized  the  motion  in  the  gasless  case  (reduced  maximum 
radius  A  *  0.620;  reduced  period  6t  *  1.135;  reduced  velocity 
p  »  0.547;  "shape  constant"  B  ■  0.625)  are,  in  the  present  case, 
functions  of  both  y  and  c.  Contour  plots  of  these  functions  are 
to  be  found  in  Figures  10  through  13.  The  reduced  maximum  radius 
A  is  always  (as  might  be  expected)  less  than  that  in  the  gaslesa 
case,  and  decreases  with  decreasing  y  and/or  C.  The  bubble 
period,  on  the  other  hand,  may  be  either  greater  or  lets  than  in 
the  gasless  case;  the  characteristic  velocity  p(  -  A/5t)  is  there¬ 
fore  lower  than  the  gasless  value  ovc  .  it  of  the  range  considered, 
but  is  slightly  larger  tor  large  values  of  y  and  of  f.  The  range 
of  variation  of  A,  6t,  and  \i  o^er  the  range  of  interest  is  consider¬ 
able;  about  17%,  9%,  and  12%,  respectively.  On  the  other  hand, 
the  "ehaoe  parameter*  B  ie  fairly  Insensitive  to  variations  in  y 
and/or  f ;  the  ovei ill  range  is  only  about  3%. 

In  order  to  evaluate  the  error  caused  by  neglecting  the 
inertial  effects  of  the  qas,  tha  same  calculations  were  then 
repeated,  out  using  the  i.cn-diaensxonal  fonts  of  equations  (VI-14) 
and  (Vi-'  ,  rather  than  (VI-18)  and  (VI-19)  ,  and  setting  (the 


45 


I 


FIGURE  8:  THE  IDEAL  GAS  MODEL  -  THE  DIMENSIONLESS  RADIUS-TIME 
RELATION  FOR  y  -  4/3  AND  RELATIVE  INTENSITY  C  •  3000 
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FIGURE  Ji 


THE  IDEAL  GAS  MODEL  -  THE  BUBBLE  ENERGY  BALANCE  rOR 
y  »  4/3  AND  RELATIVE  INTENSITY  C  -  3GC0 
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initial  gas  density)  equal  to  Pw  (the  density  of  water) .  The 
results  were  almost  identical  to  the  previous  case,  as  shown 
in  Figures  14  and  15.  Figure  14  ia  a  contour  plot  of  the  per¬ 
centage  deviaticn  in  y  between  the  two  cases  and  Figure  15  that 
in  B.  Both  y  and  B  were  consistently  lower  than  the  results 
obtained  using  pe  ■  0,  but  the  deviations  themselves  were  only 
fractions  of  a  percent.  An  examination  of  the  energy  equation 
(VI-14)  reveals,  furthermore,  that  the  maximum  bubble  radius 
is  independent  of  pe,  since  the  interface  velocity  R  is  zero 
at  the  bubble  maximum.  Thus.,  there  is  no  difference  in  A  between 
the  two  cases,  and  hence  the  deviation  in  y  (which  is  just  A/6t) 
shown  in  Figure  14  is  also  equal,  numerically,  to  the  deviation 
in  (1/6t )  . 

As  these  results  shown,  the  error  in  y  and  B  is 
almost  independent  of  y  and  depends  principally  on  c »  the 
relative  intensity  of  the  explosion.  For  explosions  at  "reason¬ 
able"  burst  depths  (that  is,  c  »  1000  or  more) ,  the  deviations 
are  less  than  one-tenth  of  one  percent;  even  considering  the 
entire  range  down  to  5  -  100,  the  doviations  are  generally  less 
than  It.  Therefore,  the  postulate  trade  earlier  that  ts  a 
inertial  effects  of  the  bubble  atmosphere  do  not  signiticantly 
affect  the  overall  motion  appears  to  be  verified. 
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VII 


REAL-STEAM  EFFECTS 


In  the  previous  section,  the  bubble  atmosphere  was 
treated  as  a  homogeneous  ideal  gas,  and  the  bubble  pulsation 
effects  were  found  to  depend  weakly  upon  y  (the  ratio  of  speci¬ 
fic  heats  for  the  ideal  gas)  and  ?  (the  explosion  intensity). 

It  wa3  also  shown  that  the  inertial  effects  of  the  bubble  atmos¬ 
phere  are  negligible  over  the  range  of  interest  of  these  para¬ 
meters.  In  reality,  of  course,  the  bubble  atmosphere  is  not 
an  ideal  gas;  therefore,  in  this  section,  the  effect  of  imposing 
a  more  realistic  equation  of  state  will  be  examined. 


Once  again,  we  will  consider  the  motion  of  a  homogen¬ 
eous  sphere  of  gas  which  is  initially  at  high  pressure  oscil¬ 
lating  in  an  incompressible  fluid.  The  governing  momentum 
equation  for  the  motion  is  therefore  the  same  as  used  previously 
(see  equations  V-12  and  VI-14): 


R 


1 

5 


H 


(VII  -1) 


This  equation  may  readily  be  integrated  numerically  if  the 
relationship  between  bubble  pressure  (Pg)  euid  bubble  volume 
is  known.  In  section  V,  it  was  assumed  that  the  bubble  pressure 
was  always  identically  xero;  in  section  VI,  tha  adiabatic 
expansion  of  an  ideal  gas  was  used  as  tha  model  for  ?3>  In 
order  to  improve  the  treatment  still  further  for  the  nuclear 
case,  the  equation  of  state  data  fcr  HjO  collected  by  Bjork, 
Kreyenhagen,  and  Wagner  (1969)  was  used. 


As  was  done  in  section  VI,  the  bubble  is  taken  as 
homogeneous,  the  motion  starts  from  rest,  and  the  expansion  is 
assumed  to  be  adiabatic  and  i sen tropic.  Ths  initial  state 
was  taken  as  that  corresponding  to  the  injection  of  the  bubble 
energy  (EQ)  into  a  spherical  volume  of  water  of  radius  equal 
to  that  of  the  "equivalent  TNT  charge”,  tilth  the  equation  of 
state  specified,  the  only  reeieining  variable  is  ths  hydrostatic 
pressure,  which  determines  the  "relative  explosion  intensity" 
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(C)  discussed  earlier.  The  equation  of  Btate  data  is  in 
graphical  and  tabular  form  and  is  therefore  rather  difficult 
to  use  for  calculations;  consequently,  only  one  case  was 
actually  computed,  corresponding  to  C  -  3270.  This  value  is 
near  the  middle  of  the  range  of  interest,  and  therefore  the 
relative  deviations  of  the  results  from  the  cases  described 
earlier  should  be  typical  of  that  to  be  expected  from  explo¬ 
sions  at  nominal  burst  depths.  The  numerical  results  are  (in 
dimensionless  form) : 

A  «*  0.6048 

<5t  =  1.1403 

p  =  0.5304 

B  »*  0.6198 

The  deviations  of  the  "gasless"  values  obtained  in  section  v 
from  the  above  are  +2.5%,  -0.51,  +0.8%,  and  +3.1%  respectively; 
that  is,  the  gasless  calculation  tends  to  slightly  overestimate 
A,  u,  and  B  and  to  underestimate  5t,  the  bubble  period.  These 
results  for  steam  may  also  be  compared  with  the  ideal  gas 
results  obtained  in  section  VI.  The  deviation  of  the  ideal  gas 
values  from  the  steam  values  ere  shown  in  Figure  16  as  functions 
of  y,  the  ideal-gas  adiabatic  exponent.  As  can  be  seen,  the 
deviations  in  all  four  quantities  are  under  1%  for  y  between 
about  1.29  and  1.38;  the  best  overall  value  would  appear  to  be 
about  4/3.  That  is,  for  this  case  at  least,  the  errors 
involved  in  replacing  the  real  steam  equation  of  state  with 
an  ideal  gas  with  y  *  4/3  are; 

A A  -  0.19% 

a6t  »  -0,14% 
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AB  »  -0 . 244 


The  conclusion  that  must  be  drawn  is  that  the  ideal 
yas  treatment  using  y  *  4/3  is  an  adequate  approximation  to  the 
real -steam  equation  of  state  as  far  as  the  bubble  dynamics  is 
concerned  As  will  be  seen  later  on,  errors  of  this  magnitude 
arc  entirely  unimportant  in  comparison  with  other  effects  and 
are,  in  fact,  of  about  the  same  order  as  those  introduced  by 
neglecting  the  inertial  effects  of  the  bubble  atmosphere  which, 
as  has  been  shown,  are  insignificant. 


viii.  fcuastt  \r^o*PHZ&*  dtssajucs 

So  fai  ,  the  bubble  interior  has  >»eer  treated  ae  if  it 
were  homogeneous  so  that  the  internal  pressure  is  independent 
of  position  within  the  buocie .  As  was  discussed  previously, 
this  necessarily  implies  that  the  internal  velocity  distribu¬ 
tion  is  given  by: 

u  -  R  (£}  for  r  <  R  (VIII-1) 

at  any  time.  In  fact,  however,  it  ,iay  easily  be  shown  that  the 
bubble  will  not  remain  homogeneous,  but  that  conditions  within 
the  bubble  will  vary  with  position  as  well  as  with  time.  In 
the  high-explosive  case,  the  bubble  may  be  roughly  homogeneous 
initially,  but  there  is  certainly  no  guarantee  that  it  will 
remain  so;  in  the  nuclear  case,  the  bubble  is  markedly  inhomo¬ 
geneous  initially.  Accordingly,  in  this  section  the  effects 
of  bubble  non-homogeneity  will  be  examined  in  an  approximate 
way  to  determine  the  consequences  as  concerns  the  overall  water 
motion.  As  before,  the  water  outside  the  bubble  is  treated  as 
incompressible,  but  the  actual  detailed  motion  within  the  bubble 
will  be  take’'  into  account  ir.  determining  its  radius- rime  his¬ 
tory. 


The  motion  of  the  gn6  within  the  bubble  comprises  a 
fully  non-linear  time-dependent  problem  in  compressible  s 
dynamics  and  is  therefore  governed  by  Euler's  equations: 
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P  e  f  ( r> ,  H  i 


(VIII-5) 


where 

u  •  velocity 
p  ■=  density 
P  =  pressure 

u2 

E  -  total  energy  per  unit  mass  =  ^  +  H 
H  =  internal  energy  per  unit  mass 
r  =  radius 
t  =  time 


Equations  Vlll-2  through  VIII-4  represent,  respectively, 
the  principles  of  mass,  momentum,  and  energy  conservation;  equation 
VI II- 5  is  the  equation  of  state  for  the  gas.  To  solve  these  equa¬ 
tions,  it  is  convenient  to  transform  them  to  a  Lagrangian  formula¬ 
tion  (that  is,  coordinates  which  move  with  the  flow,  rather  than 
remaining  fixed  in  space) .  Also,  as  was  seen  in  the  proceeding 
section,  an  ideal-gas  treatment  of  the  bubble  interior  is  adequate, 
fin a  therefore  the  equation  of  state  (V1I1-5)  is  simply: 

P  =  (y-1 ) pH 

where  y  will  be  taken  as  4/3. 

The  resulting  system  of  equations  must  be  solved  sub¬ 
ject  to  prescribed  initial  distributions  of  p,  u,  and  H,  and  to 
a  boundary  condition  at  the  bubble  perimeter  which  matches  the 
incompressible  external  motion  solved  earlier  and  which  may  be 
derived  from  the  momentum  equation  (VI-15) ,  which  defines  the 
acceleration  at  the  wall  in  terms  of  the  local  internal  pressure, 
velocity,  and  gae  density. 


To  solve  the  resulting  system  of  equations  a  finite- 
difference  method  was  employed  which  is  essentially  the  same  as 


s-iSises^ 


that  described  by  Mader  and  Ga^e  (1967).  The  hydrodynamic  code 
is  «n  explicit  forward- time  I<a*i  angian  t  lnlte-dif  f  erenc*  scheme 
employing  the  Von  Neumann- Richtm„ er  artificial  viscosity  techni¬ 
que  to  avoid  instabilities  associated  with  compression  waves. 

Thu  details  oi  the  method  have  been  extensively  described  by 
Mader  and  will  therefore  not  be  repeated  here.  Two  cases  were 
actually  calculated  using  this  method;  computer  time  limitations 
precluded  a  more  extensive  investigation. 

For  the  first  problem  run,  the  bubble  interior  was 
assumed  to  be  initially  homogeneous:  p  was  taken  as  the  same  as 
that  of  water,  end  the  "explosion  intensity"  c  wbb  3000.  The 
gaseous  region  within  the  bubble  waB  dividea  into  100  Lagrangian 
computational  cells;  as  discussed  above,  the  effects  of  the  watt 
outsidt  are  taken  into  account  by  proper  specification  of  the 
boundary  conditions  at  the  bubble  edge.  The  beginning  of  the 
bubble  expansion  caused  a  rarefaction  to  form  at  the  bubble 
edge  and  to  propagate  toward  the  interior;  the  resulting  pres¬ 
sure  Wave  continued  to  oscillate  between  the  bubble  center 
and  bubble  edge  throughout  the  rest  of  the  motion.  Several  of 
these  interior  pressure  wave  reverberations  occurred  during  the 
firct  expansion,  superimposed  on  the  general  overall  pressure 
decline.  The  Integration  was  carried  out  through  15,000  compu¬ 
tational  steps,  which  required  about  8  minutes  of  CDC-660C 
computer  time.  Calculations  were  carried  out  well  beyond  the 
point  of  maximum  expansion  and  into  the  beginning  of  the 
bubble  collapse.  Although  the  governing  equations  no  longer 
imply  symmetry  of  the  radius-time  curve,  the  computed  results 
themselves  are  symmetric  to  within  small  fractions  of  a  percent. 
Accordingly,  the  bubble  period  was  taken  as  juot  twice  the 
time  to  maximum  expansion. 

The  maximum  bubble  radius  in  dimensionless  form 
computed  in  this  way  was  slightly  less  than  that  obtained 
assuming  a  homogeneous  bubble  atmosphere,  since  the  kinetic 
energy  oi.  the  bubble  atmosphere  did  not  fall  to  zero  at  the 
maximum;  the  oscillating  pressure  disturbance  mentioned  above 
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MODEL 

Dynamic 

Ideal  Gas 

Ideal  Gas 

Bubble 

With 

Without 

Gaslese 

REVOLT 

Interior 

Inertial  Effects  Inertial  Fffectc 

Model 

A 

0.6053 

0.6C56 

0.6066 

0.620 

6x 

1.1392 

1.1394 

1.1389 

1.135 

U 

0.5314 

0.5315 

0.5317 

0.547 

2 

0.6177 

0.6178 

0.6180 

0.625 

As  can  be  seen,  the  consideration  cf  the  interior  gas  dynamics 
made  virtually  no  difference  to  the  results,  and  consequently 
the  homogeneous  approximation  appears  to  be  quite  adequate  to 
describe  the  notion.  Although,  in  the  dynamic  calculation, 
fairly  large  pressure  fluctuations  around  the  mean  value 
occurred  at  the  bubble  boundary,  a  sufficient  number  of  fluc¬ 
tuations  occurred  that  the  average  overall  effect  was  about 
the  same  as  if  there  had  been  no  such  fluctuations. 

In  the  caee  juet  diBcussed,  however,  the  bubble  was 
assumed  to  be  initially  homogeneous.  As  was  discuseed  earlier 
a  nuclear  explosion  would  be  expected  to  produce  an  "initial 
bubble"  that  is  non-homogeneous ,  that  is,  hotter  and  less 
dense  at  the  center  than  at  the  edge.  To  evaluate  the  effect 
of  initial  non-homogeneity  of  the  bubble,  the  same  Lagrangian 
hydrodynamic  code  was  used  as  in  the  previous  case,  and  all 
input  conditions  were  the  same,  except  that  for  the  innermost 
fourth  of  the  initial  bubble  radius,  the  density  was  set  at 
one- tenth  of  the  previous  value  and  the  internal  energy  per 
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unit  maes  (which  is  proportional  to  temperature!  was  multiplied 
by  tan.  For  the  outer  three-fourths,  the  initial  conditions 
were  tie  Bane  as  in  the  previous  case. 

The  second  problem  was  not  carried  out  in  time  as  far 
as  the  first  bubble  maximum,  since  the  resultB  even  at  early 
times  were  essentially  identical  (as  far  as  the  bubble  radius- 
time  relation  was  concerned)  with  those  in  the  previous  rase. 
The  deviation  never  exceeded  0,04%  in  magnitude,  and  declined 
slowly  toward  zero  as  the  computation  went  on.  The  computa¬ 
tion  was  stopped  at  t  ■  0.17,  at  which  time  the  bubble  radius 
was  75%  of  the  maximum  value,-  at  that  point,  the  deviation 
in  bubble  radius  between  the  two  cases  had  dropped  to  C.C27i. 
The  overall  conclusion  which  must  necessarily  be  reached  from 
both  of  these  calculations  is  that  the  explicit  treatment  of 
the  bubble  atmosphere  dynamics  does  not  significantly  alter 
the  results  for  the  bubble  radius-time  relation  from  those 
obtained  using  the  homogeneous  model  discussed  in  section  VI, 
and  that  this  is  true  regardless  of  the  initial  energy  distri¬ 
bution  within  the  bubble. 
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tx.  COMPARISONS  WITH  OTHER  RESULTS 

To  evaluate  the  adequacy  of  the  models  discussed  in 
the  preceedinq  sections,  it  is  nscesBary  to  consider  experi¬ 
mental  results.  As  was  just  shown,  the  effects  upon  the  radius- 
time  relation  of  a  non-homogeneous  dynamic  bubble  interior  are 
extremely  alight.  Furthermore,  over  the  range  of  interest,  the 
effects  of  the  deviations  of  real  steam  behavior  from  y  -  4/3 
ideal  gas  behavior  are  unimportant,  and  the  inertial  effects  of 
the  bubble  atmosphere  are  likewise  negligible.  Of  the  various 
models  considered,  only  the  results  of  the  most  primitive  (that 
is,  the  "gasless"  model  of  section  V}  differ  to  any  significant 
extent  from  the  rest,  and  even  that  case  is  generally  within  a 
few  percent  of  the  others.  The  general  conclusion  that  must  be 
drawn  is  that  the  bubble  motion  is  quite  insensitive  to  the 
behavior  of  the  interior;  that  is,  that  the  motion  is  dominated 
by  the  water  flow  outside  the  bubble  rather  than  that  of  the 
gas  inside. 

In  all  cases  considered,  the  water  outside  the  bubble 
was  assumed  to  be  incompressible.  This  assumption  is  appro¬ 
priate  in  view  of  the  overall  purpose  of  thiB  investigation  and 
has  further  been  shown  by  various  plausibility  arguments  to  be 
at  leaBt  approximately  correct  except  for  short  time  periods  very 
near  the  bubble  minimum.  None  the  less,  it  is  certainly  worth¬ 
while  to  determine  the  extent  to  which  the  ideal  gas  -  incompres¬ 
sible  water  model  represents  observed  explosion  bubble  behavior. 

An  enormous  amount  of  experimental  data  has  been 
acquired  over  the  years  concerning  underwater  explosions.  Even 
so,  precise  measurements  of  the  bubble  motion  are  extremely 
difficult  to  obtain.  For  large  explosions  in  the  field,  the 
only  bubble  parameter  which  may  be  readily  determined  is  the 
period  of  oscillation,  that  is,  the  time  interval  between 
the  emission  of  the  primary  shockwave  and  the  "bubble  pulse" 
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piessura  wave  emitted  jc  the  first  reoompress ion  For  high 
explosives  fired  deep  enough  that  the  facsonce  ui  tho  free 
surface  does  not  significantly  affect  the  motion,  the  data  may 
b«j  correlated  by  an  empirical  expression  of  the  form  (Cole , 
1948)  : 


T (sec) 


WT  (lh.  ) 

TTTt~ 


. ) 
Tj 


1/3 


(1X-1) 


where  WT  is  the  charge  weight  in  pounds,  Z  is  the  total  hydro¬ 
static  pressure  measured  in  feet  of  water  (that  is,  the  depth 
of  burst  plus  34  feet  representing  air  pressure) ,  and  X  depends 
upon  the  .  .rticular  type  of  explosive  used.  For  TNT,  for 
example , 


T(sec)  *  4.36 


WT (lb. ) 

TT7E7T 


1/3 


(TX-2 ) 


will  predict  observed  bubble  periods  within  a  few  percent  over 
a  wide  range  of  charge  weights  and  burst  depths.  If  the  initial 
bubble  energy  for  TNT  is  taken  as  47%  of  the  total  explosion 
energy  (that  is,  the  total  minus  the  observed  53%  shockwave 
fraction) ,  the  total  energy  per  unit  explosive  mass 

is  assumed  to  be  4.2  *  106  joulea/kilograms  (see  section  II), 
and  water  density  la  taken  as  1000  kg/m3,  equation  (IX-2)  may 
be  expressed  in  dimensionless  form: 
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(IX- 3) 


that  is, 


6  t 


1.126 


which,  as  was  shown  in  section  VI  (see  Figure  11),  agrees  quite 
well  with  the  vaiues  obtained  from  the  homogeneous  ideal-gas 
model.  Over  the  region  of  interest,  the  ideal  gas  results  for 
y  *  1.25  (the  valuo  suggested  by  Snay  (1957)  as  being  most 
appropriate  for  TNT)  are  illustrated  in  Figure  17 
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FIGURE  17:  COMPARISON  OF  THE  "AVERAGE"  EXPERIMENTAL 
DIMENSIONLESS  OSCILLATION  PERIOD  FOR  TNT 
WITH  IDEAL  GAS  MODEL  RESULTS 
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over  the  range  £  ■»  1000  to  C  ■  30,000  (the  region  in  which  rooet 
oi  the  experimental  data  was  acquired)  la  always  less  than  0.15%, 
far  less  than  the  scatter  in  the  experimental  data  itself.  This 
rather  remarkable  agreement  between  the  ideal-gas  results  and 
experimental  measurements  for  TNT  is  probably  somewhat  fortui¬ 
tous,  but  certainly  indicates  that  the  model  can  adequately  re¬ 
present  experimental  .esults  for  tha  bubble  oscillation  period. 

For  nuclear  explosions,  experimental  data  is  extremely 

sparse.  The  only  available  direct  measurement  concerning  nuclear 

explosion  bubble  motion  is  the  bubble  period  measured  hydro- 

acoustically  at  the  1955  WIGWAM  test.  Operation  WIGWAM  consisted 

of  the  detonation  of  a  30  kiloton  nuclear  device  (total  energy 
14 

release  1.2  *  10  joules)  at  a  depth  of  610  meters  in  very  deep 
water  in  the  open  sea.  The  first  bubble  period  was  2.88 
seconds; if  the  fluid  density  is  taken  as  1025  kg/m^  (an  "average" 
value  for  seawater)  and  the  initial  bubble  energy  is  taken  as 
40%  of  the  total  yield  as  discussed  in  section  II,  the  dimension¬ 
less  bubble  period  is: 

6 t  «  1.137 

which  may  be  compared  with  the  results  of  the  y  *■  4/3  ideal-gas 
model  (see  Figure  18).  The  "reduced  explosion  intensity"  £ 
for  WIGWAM  was  taken  as  500,  based  upon  the  assumptioh  of  an 
initial  bubble  energy  density  equal  t<  that  of  TNT.  The  devia¬ 
tion  between  the  two  is  about  1%,  which  is  within  the  precision 
of  the  WIGWAM  oscillation  period  measurement. 

There  are,  of  course,  no  measurements  of  maximum  bubble 
radii  for  underwater  nuclear  explosions.  Even  for  high  explo¬ 
sive*,  far  less  data  exists  for  the  bubble  radius  than  for  the 
bubble  period.  What  data  is  available  may,  however,  be  corre¬ 
lated  by  the  empirical  formula: 
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which  is  analogous  to  the  empirical  equation  for  the  bubble 
period  (1X-1) .  For  TNT,  J  ■  12.6  fits  most  of  the  available 
data  within  a  few  percent.  This  observed  result  may  be  trans¬ 
formed  to  non-dimensional  form  in  the  same  way  as  was  done  with 
the  TNT  bubble  period?  the  result  is: 


A  *  0.58 


which  may  be  compared  with  the  ideal-gas  model  results  using 
y  =  1.25  (see  Figure  19).  Again,  the  agreement  is  fairly  good. 

Many  of  the  available  meacur;;.....-  ui  bubble  motion 
were  taken  on  a  laboratory  scale,  using  very  tiny  charges  in 
special  test  chambers  in  which  the  air  pressure  may  be  reduced 
to  provide  Froude  scaling,  that  is,  similitude  of  bubble  buoy¬ 
ancy  effects  isee,  for  example,  Snay ,  1964;  Pritchett ,  1966) . 

The  purpose  of  these  tests  was,  in  general,  to  simulate  the 
buoyant  bubble  migration  characteristics  of  large  explosions. 

For  the  present  purpose,  however,  much  of  this  data  must  be 
regarded  with  some  suspicion.  The  reason  is  that,  at  low 
ambient  hydrostatic  pressures,  the  local  pressure  in  the  water 
near  the  bubble  interface  may  drop  below  vapor  pressure  near 
bubble  maxima  and  the  surface  of  the  bubble  will  then  boil. 

Figure  20  shows  the  pressure  distribution  in  the  water 
for  an  explosion  at  (  «  3000  and  y  ■  4/3.  The  ideal-gas  model 
of  section  VI  was  used  for  this  calculation,  and  the  inertial 
effects  of  the  gas  were  neglected.  The  pressure  distribution  was 
then  obtained  using  equation  V-12.  As  can  be  seen,  at  early 
times,  the  pressure  declines  with  increasing  distance,  but  once 
the  internal  pressure  drops  below  PH  (the  hydrostatic  pressure) , 
the  point  of  lowest  instantaneous  pressure  is  directly  adjacent 
to  the  bubble  interface.  Thus,  if  the  interior  pressure  at  the 
maximum  predicted  by  this  theory  is  less  than  the  water  vapor 
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FIGURE  19:  COMPARISON  OF  THE  "AVERAGE"  EXPERIMENTAL 

DIMENSIONLESS  MAXIMUM  BUBBLE  RADIUS  FOR  TNT 
WITH  IOEAL  GAS  MODEL  RESULTS 
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pressure  (which  depends  on  temperature) ,  the  water  at  the  bubble 
boundary  will  begin  to  boiJ ,  increasing  the  amount  of  gas  in 
the  bubble,  ite  maximum  size,  and  its  period  of  oscillation.  For 
explosions  at  normal  hydrostatic  pressures  (greater  than  one 
atmosphere) ,  of  course,  this  effect  does  not  occur.  If  the 
water  is  cooled  sufficiently,  bubble  boiling  will  not  occur  even 
for  hydrostatic  pressures  less  than  one  atmosphere.  The  hydro¬ 
static  pressure  at  which  bubble  boiling  at  the  maximum  would  be 
expected  to  begin  as  a  function  of  water  temperature  is  shown  in 
Figure  21.  This  result  should,  of  course,  be  taken  as  approximate; 
it  was  derived  using  tins  >  =  4/3  ideal  gas  bubble  model  assuming 
an  initial  bubble  energy  density  equal  to  that  of  TNT.  Actually, 
Figure  21  is  quite  conservative;  experimental  results  indicate 
that  significant  effects  upon  bubble  behavior  do  not  occur 
unless  hydrostatic  pressures  are  about  a  faccor  of  four  or  so 
smalier  than  indicated  (Snay ,  1964)  .  The  reason  is  simply  that, 
although  the  bubble  interior  pressure  may  drop  slightly  below 
vapor  pressure  ai  the  bubble  maximum,  the  amount  of  water  actually 
vaporized  is  fair  y  small  due  to  the  >  sly  short  bubble 

oscillation  period,  unless  the  interior  pressure  is  well  below 
vapor  pressure  foi  a  substantial  part  of  the  oscillation  cycle. 

A  number  of  experiments  using  the  "vacuum  tank  techni¬ 
que"  were  performed  at  the  Naval  Radiological  Defense  Laboratory 
during  the  early  ISaO’s  using  a  very  high  energy  density  elec¬ 
trically  exploded  metallic  wire  as  the  energy  source  (Buntzen, 

1964;  Pritchett,  1966).  Such  an  explosion  produces  no  permanent 
gases,  and  therefore  the  steam  bubble  generated  is  similar  to 
that  of  an  unde; water  nuclear  explosion  (Buntzen,  1961).  Some 
of  theae  tests  occurred  in  the  "danger  zone"  in  which  the  bubble 
boiling  might  be  expected;  of  those  which  were  not,  however,  the 
bubble  radii  and  oscillation  periods,  as  measured  by  high-speed 
photography,  agree  fairly  well  with  results  predicted  using  the 
ideal-gas  model.  The  scatter  in  this  data  was,  however,  con¬ 
siderable;  due  to  the  small  size  of  the  explosions,  bubble  periods 
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(atmospheres) 


were  of  order  a  few  tens  of  milliseconds  and  bubble  radii  were 
typically  only  a  few  centimeters.  Furthermore,  thu  tot.al 
explosion  energy  was  difficult  to  reproduce,  and  could  only  be 
determined  within  10%  or  so.  None  the  less,  these  results  also 
tend  to  support  the  ideal-gas  model  for  bubble  pulsation. 

Other  measurements  than  the  maximum  radius  and  bubble 
period  have,  of  cource,  been  made.  The  scatter  in  the  measured 
results  is  considerable,  however,  due  to  the  great  difficulty 
in  experimentally  determining  the  bubble  radius-time  relation. 

The  radius-time  curve  computed  using  the  ideal-gas  model  falls 
well  within  the  band  of  experimental  scatter.  Measurement  of 
subtle  parameters  such  as  B(the  "shape  factor")  would  of  cv.  *ue 
be  extremely  difficult,  and  has  never  been  attempted,  to  the 
author's  knowledge.  The  "characteristic  velocity"  p,  on  the 
other  hand,  may  be  determined  if  the  bubble  maximum  radius  and 
oscillation  period  are  known.  The  "average"  value  which  may  be 
derived  from  the  empirical  TNT  equations  (IX-2  and  IX-4)  is  0.51, 
which  falls  in  the  middle  of  the  range  of  values  derived  from  the 
ideal  gas  model.  As  has  been  pointed  out  by  Snay  (I960)  and 
others,  although  the  empirical  coefficients  J  and  K  depend 
upon  explojive  type,  the  ratio  j/K  *.1  relatively  constant  at 
about  2.9  within  a  few  percent,  which  implies  p  ■  0.51. 

As  has  been  seen,  the  available  experimental  data,  while 
generally  verifying  the  present  results  which  assume  incompres¬ 
sible  water,  are  sufficiently  scattered  that  precise  comparisons 
are  simply  not  possible.  In  recent  years,  a  few  theoretical 
calculations  of  one-dimensional  bubble  motion  which  include  the 
effect  of  water  compressibility  have  appeared  in  the  literature. 
Two  of  these  (Kot ,  1964  and  Bjork,  Kreyenhagen  and  Wagner,  1969) 
involved  nuclear  explosions,  and  two  others,  (Phillips  and  8nay, 
1968  and  Mader ,  1971)  represented  smaller  conventional  explosions. 
Kot  considered  three  cases,  all  one-kiloton  nuclear  explosions 
at  burst  depths  of  91,  305,  and  610  meters.  The  equation  of 
state  data  used  war,  regrettably,  rather  primitive.  Bjork,  et. 
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al .  considered  only  one  case  —  a  10  kiloton  nuclear  explosion 
at  1830  meters.  Unfortunately,  the  one -dimensional  calculation 
was  terminated  at  6.9  microseconds ;  a  relatively  poorly-resolved 
two-dimensional  method  was  used  thereafter,  but  even  thiB  cal¬ 
culation  ended  at  140  milliseconds,  so  that  the  bubble  never 
reached  its  first  maximum.  Phillips  and  Snay  report  the  results 
of  five  calculations  of  one-pound  TNT  explosions  at  burst  depths 
oi  24,  152,  305,  4220  and  6822  meters.  Mader ' s  results  repre¬ 
sent  the  explosions  of  half-pound  spherical  Tetryl  charges  (a 
conventional  high  explosive  similar  in  characteristics  to  TNT) 
at  hydrostatic  pressures  corresponding  to  burst  depths  of  89, 

7  35,  4‘ilO  and  46,600  meters.  In  reality,  of  course,  the  average 
depth  of  the  ocean  is  only  about  4000  meters  and  the  deepest 
point  is  about  11,000  meters  (at  the  bottom  of  the  Mindanao 
Trench)  bo  Mader ' s  deepest  case  must  be  regarded  as  somewhat 
hypor.hetical .  The  range  of  depths  of  interest  for  underwater 
explosions  generally  does  not  extend  below  1000  meters  or  so. 

The  maximum  bubble  radii  reported  by  Kot,  Phillips 
and  Snay,  and  Mader  (excluding  his  deepest  case  which  is  off- 
scale)  are  shown  in  dimensionless  form  as  functions  of  the 
relative  explosio:  intensity  c  in  Figure  22.  Also  shown  are 
the  results  of  the  ideal  gas  model  neglecting  gas  inertial 
effects  for  y  ranging  from  1,15  to  1.40.  The  scatter  in  the 
"band"  of  compressible  -water  computational  results  1b  about 
10%,  but,  as  can  be  seen,  the  overall  agreement  is  fairly  good. 
Similarly,  the  results  of  Kot,  Phillips  and  Snay,  and  Mader 
for  the  "characteristic  velocity"  u(  ■  A/6t)  are  compared  with 
the  ideal-gas-model  results  in  Figure  23.  Also  included  is 
the  result  of  an  earlier  calculation  by  Keller  and  Kolodner 
(1953)  in  which  the  water  was  considered  compressible,  but  the 
bubble  interior  was  treated  as  a  homogeneous  ideal  gas  with 
y  *  1.25.  Once  again,  the  agreement  of  the  data  with  the 
incompressible  results  appear  as  good  as  the  agreement  between 
one  set  cf  data  and  another. 


COMPARISON  OF  THE  DIMENSIONLESS  MAXIMUM  BUBBLE  RADII  (/)  PREDICTED  BY 
COMPRESSIBLE-WATER  CALCULATIONS  WITH  IDEAL  GAS  -  INCOMPRESSIBLE  WATEF 
RESULTS 


RESULTS 


Although  both  the  experimental  measurements  discussed 
earlier  and  the  theoretical  results  which  include  the  effect 
of  water  compressibility  tend  to  agree  with  the  incompressible- 
water  calculations,  the  scatter  in  both  is  sufficiently  great 
that  no  clear  choice  based  upon  these  results  alone  can  be  r ade 
among  the  various  models  discussed  in  the  previous  sections.  The 
ideal-gas  model  results  agree  quite  well  wi^h  experiments  and  with 
the  compressible-water  calculations  and,  as  has  been  shown, 
further  elaboration  of  that  model  (to  include  gas  inertial 
effects,  real-steam  behavior,  and/or  bubble  inhomogen'Uty) 
alters  the  results  obtained  by  only  fractions  of  a  peicer.t. 

The  d£  'iations  of.  the  ideal-gas  incompressible-water  model 
results  from  compressible-water  calculations  and  from  experi¬ 
mental  measurements  are  at  least  as  great  as  the  "improvement" 
that  could  in  principle  be  obtained  by  using  a  more  elaborate 
model  for  the  bubble  interior. 

It  is  therefore  recommended  that  the  ideal  gas  model 
neglecting  gas  inertial  effects  be  adopted  for  the  purpose  of 
making  incompressxble-water  calculations  of  bubble  motion.  The 
additional  complications  introduced  by  using  a  more  complex 
model  are  simply  not  warranted  in  view  of  the  very  slight  gain 
in  accuracy,  particularly  since  the  assumption  of  incompressible 
water  itself  introduces  errors  of  at  least  comparable  size.  As 
was  seen  in  section  VII,  fox  steeun  bubbles  produced  by  nuclear 
explosions,  y  *  4/3  seems  to  agree  best  with  the  real-ttaam 
results?  for  h-gh-explosive  calculations,  a  somevnat  lower  value 
(such  as  1.25)  should  probably  be  used. 
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APPENDIX  A:  SYMBOLS  AND  ABBREVIATIONS 


A 

B 


C 


9 

H 

I 

I* 


J 

K 

M 

P 

P. 


instantaneous  bubble  surface  area 

"shape  factor"  for  bubble  oscillation  cycle  ■ 


speed  of  sound  in  water 

detonation  wave  speed  within  the  explosive 
total  energy  per  unit  mass  -  ^  u2  +  H 
bubble  oscillation  energy 
internal  energy  of  bubble  atmosphere 
initial  oub^le  internal  energy 
kinetic  energy 

kinetic  energy  of  bubble  atmosphere 

kinetic  energy  of  the  water  outside  the  bubble 

buoyant  force 

acceleration  of  gravity 

internal  energy  per  unit  mass 

buoyant  momentum 

dimensionless  buoyant  momentum 

empirical  "radius  coefficient" 

empirical  "period  coefficient” 

Mach  number 
pressure 

initial  bubble  pressure 
bubble  pressure 
hydrostatic  pressure 

chemical  energy  released  per  unit  explosive  mass 
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r 


radial  coordinate 


i 


z 

6  T 

y 

x 

x 

x 

o 

A 

W 


instantaneous  bubble  radius 

instantaneous  bubble  interface  velocity 

instantaneous  bubble  interface  acceleration 

initial  bubble  radius 

maximum  bubble  radius 

time 

reaction  time  of  explosive  charge 
bubble  oscillation  period 
velocity 

characteristic  velocity  for  bubble  oscillation  ■ 

■WT 

bubble  volume 

initial  babble  volume 

work  done  against  external  forces 

charge  weight  measured  in  pounds 

yield;  total  explosion  energy 

axial  coordinate  variable 

hydrostatic  pressure  (PH)  measured  in  feet  cf  water 
dimensionless  bubble  oscillation  period 
ideal-gas  adiabatic  exponent  (ratio  of  specific  heats) 
dimensionless  instantaneous  bubble  radius 
dimensionless  instantaneous  bubble  interface  velocity 
dimensionless  initial  bubble  radius 
dimensionless  maximum  bubble  radius 

dimensionless  bubble  oscillation  characteristic  velocity 

(A/6  t ) 
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it  3.14153255... 

p  density 

PF  explosive  packing  density;  initial  bubble  density 

Pg  instantaneous  bubble  interior  density 

Pw  water  density 

t  dimensionless  time 

t  dimensionless  explosion  "intensity”. 
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